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ABSTRACT 

We examine the spin-up of low luminosity, low mass X-ray binaries (LMXBs) to 
millisecond pulsars (MSPs). In the conventional spin- up model of the Ghosh & Lamb 
type, where the stellar magnetic field interacts with the Keplerian accretion disk, MSPs 
could be produced from LMXBs if magnetic field B* < 10 8 (M/10 16 5/s) 1 / 2 G, where M 
is the mass accretion rate. However, for M < M c ~ 10 16 g/s accretion is likely to occur 
via a quasi-spherical flow with a sub-Keplerian rotation. The sub-Keplerian rotation 
rate is smaller than the Keplerian rate by a factor ~ 2 — 10. As a consequence, the 
spin-up of LMXBs produces pulsars with spin periods longer by a factor ~ 2 — 10 than 
those with a Keplerian accretion disk. The observed MSPs could be produced only 
for B* < 10 G even when M ~ M c ~ W ie g/s. This suggests that the low luminosity 
LMXBs with M < M c would not be able to spin-up to the observed MSPs. This rules 
out any undetected populations of persistent, low luminosity LMXBs and potentially a 
significant fraction of the atoll sources as a possible source population of the observed 
MSPs. If a large number of undetected, persistent, low luminosity LMXBs do exist, 
they could produce MSPs near the pulsar death line with intrinsic electromagnetic 
luminosity ^ 10 30 erg/s. The observed MSPs could possibly arise from a population of 
soft x-ray transients (SXTs) containing neutron stars although this is not supported 
by current estimates of the outburst rate or numbers of such systems. Accretion 
induced collapse of low magnetic field white dwarfs remains a possible channel for MSP 
formation. 



Subject headings: accretion, accretion disks — pulsars: general — stars: magnetic fields 
— X-rays: stars 



1. Introduction 

The origin of millisecond pulsars (MSPs) remains an outstanding problem (e.g. Bhattacharya 
& van den Heuvel 1991). A number of channels of MSP formation have been discussed. In the 
standard model, magnetized, low mass X-ray binaries (LMXBs) spin-up to MSPs through mass 
accretion (e.g. Frank et al. 1992). Several statistical analyses, however, have questioned whether 
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the known LMXB populations could produce the observed MSPs (Grindlay & Bailyn 1988, 
Kulkarni & Narayan 1988, Bailyn & Grindlay 1990, Kulkarni 1995). 

The problem of MSP birth rate is briefly summarized as follows (e.g. Grindlay 1995). 
The MSP lifetime due to the electromagnetic dipole emission (e.g. Shapiro & Teukolsky 1983) 
is tmsp ~ 10 10 (-B*/5 x 10 8 G) _2 (P*/5ms) 2 yr where B* is the pulsar magnetic field strength 
and P* is the spin period. For the estimated total number of MSPs (Bailes & Lorimer 1995), 
Nmsp ~ 10 5 , the observed MSPs require a birth rate of ~ Nmsp/tmsp ~ 10~ 5 on the MSP 
birth rate, ~ 3 x 1CT 6 yr~ x above the luminosity limit of 1 mJy kpc 2 . The uncertainty in this 
estimate could be large due to uncertainties in beaming of pulsar emission and the Galactic 
scale height of the pulsar distribution. For a high accretion rate M and a relatively low B*, 
the accretion flow extends close to the neutron star surface. Then, the spin-up of a slowly 
rotating neutron star to the spin period of ~ 5ms, would require mass accretion of at 
least > AM ~ /*Q*/( GM *#*) 1/2 ~ O.1M , where /* = 10 45 gem 2 is the moment of inertia, 
M* = 1.4M is the mass, = 10 6 cm is the radius of the neutron star, and 17* = 2ir/P*, 
respectively (cf. Bhattacharya & van den Heuvel 1991). Luminous Z sources have luminosities 
Lz ~ 3 x 10 37 — 10 38 erg/s close to the Eddington luminosity L^dd ~ 2 x 10 3S erg/s and their 
magnetic fields are estimated to be in the range ~ 10 9 ~ 10 G. The typical spin-up time scale for Z 
sources tz <^ GM^AM/ LzR* ~ 3 x 10 7 yr. For the observed number of Z sources Nz = 6, the 
MSP birth rate from Z sources is expected to be ^ Nz/t z ~ 2 x 10~ 7 yr~ 1 . For less luminous, 
bursting atoll sources with luminosities L ato ii ~ 10 36 — 10 37 erg/s, the typical spin-up time scale 
Tatoii ~ GM^AM I L at oiiR* ~ 10 9 yr. For the total number of atoll sources N ato ii ~ 10 2 , the 
estimated birth rate is ^ N a t u /r a toii ~ 10~ 7 yr _1 , which is within a factor ~ 2 of the Z source 
rate. These rates fall short of the MSP formation rate roughly by at least an order of magnitude 
(cf. Kulkarni 1995). 

Recent observations with RXTE (cf. summary in White and Zhang 1997) have yielded the 
long-sought evidence that LMXBs (including both Z sources and burst sources) are spun-up to 
~few ms spin periods as indicated by their kHz QPOs and the direct detection (in a few sources) 
of what appear to be underlying spin periods. Some 10 sources now appear to have ~ 300Hz 
pulsation frequencies which certainly suggests that these sources have indeed spun-up and could 
thus evolve to MSPs. Recently, the first direct detection of a true LMXB-MSP (SAXJ1808.4-3658) 
with a spin-period of 2.49ms has been reported (Wijnands and van der Klis 1998) along with the 
value for its binary period (~2h), probable companion mass (^ 0.1M Q ) and thus likely evolution 
(Chakrabarty and Morgan 1998). However, even if all the known LMXBs are found to have similar 
pulsation periods, the birth rate problem summarized above still remains to be resolved. 

The MSP formation problem could be solved if a generally unseen population of LMXBs 
with low luminosities spin-up to MSPs. Recent observations have measured low quiescent 
luminosities from recurrent transients as well as some usually bright persistent sources (Verbunt 
et al. 1994). The required number iVj of the low luminosity LMXBs with the luminosity L\ is 
Ni ~ 10 N ato u(L ato u j Li) if these sources are the solution to the MSP formation problem. Since 
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the spin-up time scale for the low luminosity population 77 ~ T a t ii{L a t ii / L{) , L\ must be at 
least ~ 0.1L ato u in order to avoid an excessively long spin- up time scale. Ni ~ 100N ato u is 
required with L; ~ 0.1L ato u. The low luminosities are either a constant for persistent sources or a 
time-averaged value for recurrent transient sources. It is crucial whether low luminosity sources 
with L\ ~ O.lLatoii, persistent or recurrent, are able to spin up to MSPs. 



2. Spin-Up of LMXBs: Standard Model with Low M 

The spin period of the neutron star evolves according to P* = — iVP 2 //*, where N is the 
torque exerted on the star by the accretion flow. In the conventional magnetized, Keplerian 
accretion disk mode (e.g. Wang 1995, Yi et al. 1997, and references therein), the accretion flow is 
truncated by the neutron star magnetosphere at a radius R Q determined by 

(Ro/Rc) 7/2 [l - {R /Rcf /2 \ ~* = 2N C /M(GM*R C ) 1 / 2 (2-1) 

where R c = (GM^P 2 /47T 2 ) 1 / 3 is the Keplerian corotation radius, N c = (~//a)B 2 R®R~ 3 , a is 
the usual viscosity parameter (Frank et al. 1992), and 7 is the parameter of order unity which 
determines the pitch of the magnetic field in the accretion disk (e.g. Wang 1995). For numerical 
values, we adopt a = 0.3 and 7 = 1 unless noted otherwise. The torque exerted on the star 

TV = (7N /6) [l - (8/7)(P /P c ) 3 / 2 ] / [l - (R /R c )^ 2 } (2-2) 

where iV = ^(GNURo) 1 ! 2 . The spin equilibrium N = occurs when R /R c = x eq = (7/8) 2 / 3 . 
This adopted magnetic torque model differs little from other phenomenological models for the 
purpose of the present discussion (Wang 1995). 

In the spin-up regime TV > or R /R c < x eq , we get from eq. (2-2) 

N w 1.3(GM*R 2 f/ 7 '(7/a) 1 / 7 Bl n M 6 / 7 . (2-3) 

In this regime, the spin period evolves simply as oc The spin-up line is 

determined by the equilibrium spin condition N = 0, 

P* « [bms^aZf 7 < 8 7 Mr 6 3/7 (2-4) 

where a_i = a/0.1, P* i8 = P*/10 8 G, and M 16 = M/10 16 g/s. Spin periods of P* < 5ms would 
be reached for B t < [3 x 10 8 G](a/ 7 ) 1 /2 m|/ 2 . 

The spin-up line corresponding to M = W ls g/s, 
which is close to the Eddington rate, is P* ~ [0.4ms]P^g 7 , which obviously implies that MSPs 
with P* ~ 5ms could easily form from high luminosity LMXBs such as Z sources with B* as 
high as ^ 3 x 10 9 C For M = 10 16 <7/s, which is relevant for atoll sources, the spin-up line is 
P* [3ms] 5^ g 7 . P* ~ 5ms requires B»^2x 10 8 G. Given the dependence of P* oc (a/7) 1 / 2 , it is 
gratifying that the uncertainties in a and 7 are not critical unless a < 0.1 and 7>1, which is 
highly unlikely (e.g. Aly & Kuijpers 1990, Wang 1995). 
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In Figure 1, we plot the spin-up lines in the vs. plane. The x's mark the positions 
of the observed pulsars with their spin periods and estimated magnetic fields adopted from 
Taylor et al. (1995). The magnetic fields have been estimated based on the assumption that 
their power is derived from the electromagnetic dipole emission at a rate E = Airl^P^/P 3 ; i.e. 
B* = (3/*c 3 P*P /87T 2 -R*) 1 / 2 . In Figure la, we show spin-up lines assuming that the accretion flow 
remains Keplerian for all M's. This plot shows that MSPs form only from LMXBs with weak 
magnetic fields < 10 8 G if M < 6 x 10 15 g/s or LMXBs with B* < 10 7 G if M ~ 10 13 g/s. A 
system with a lower M needs proportionally more time to spin-up to a MSP than does a high 
M system. The final spin period becomes the equilibrium period, which is fixed by both M 
and B*, only if the system arrives at the equilibrium period on timescales for the binary not to 
have evolved significantly. Otherwise, M would have changed significantly, in which case a new 
equilibrium spin period would be sought. For M ~ W 13 g/s, the spin-up time scale is prohibitively 
long 3> 10 10 yr and spin-up to a MSP in a time-scale less than the age of the Galaxy is impossible. 
Therefore, in the standard model, the MSP formation rate problem could be solved only if there 
are ^ 5x 10 2 N a t o u ~ 5 x 10 4 persistent, low luminosity LMXBs with M ~ 10 15 g/s (i.e. luminosity 
L~5x 1CP 2 L ato u) and magnetic fields ~ 5 x 10 7 G (see below for recurrent sources). 



3. Spin-Up by Sub-Keplerian Flows at Low Mass Accretion Rates 

The conclusion we have derived above may get dramatically worse due to the possible 
transition in accretion flows (Grindlay 1995). This is because at low M, the accretion flow is likely 
to take the form of a quasi-spherical, sub-Keplerian rotation (e.g. Narayan & Yi 1995). More 
specifically, when M falls below a critical accretion rate M c , the viscously dissipated energy in 
the accretion flow is not efficiently radiated but kept as internal heat energy of the accreted gas. 
As a result, the accretion flow is heated essentially to the virial temperature and the accretion 
flow thickens due to the internal pressure support. The internal pressure support also causes the 
rotation of the quasi-spherical flow to deviate from Keplerian to sub-Keplerian. 

The sub-Keplerian rotation of the accretion flow can be conveniently modeled as f2 = uJoQk 
where u> a < 1 where = (GM*/ 'E?) 1 / 2 is the Keplerian rotation (Narayan & Yi 1995, Yi et 
al. 1997). The sub-Keplerian rotation parameter uo a is determined essentially by the ratio of the 
internal pressures (3 = P gas /(Pgas + Pmag) where P gas is the gas pressure and P ma g is the magnetic 
pressure from the tangled, isotropic magnetic field in equipartition with the thermal energy of the 
accreted plasma; i.e. u a = [(10 — 67 s )/(97 s — 5)] 1 / 2 where 7 S = (8 — 3/3)/(6 — 3/3) is the ratio of 
specific heats in the assumed equipartition plasma (Esin 1997, cf. Narayan & Yi 1995). A large 
(3 ~ 1 (i.e. small magnetic pressure) would give u> a 1 whereas uj a ~ 0.4 for equipartition (3 ~ 0.5. 
The sub-Keplerian rotation parameter oj a weakly depends on a as long as a = 0.01 — 0.3. Recently, 
Yi et al. (1997) have suggested that the torque reversals seen in some X-ray pulsars could be 
precisely due to this accretion flow transition, which would also confirm that indeed sub-Keplerian 
flows form probably for M below a critical rate M c for accretion onto NSs since otherwise these 
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advection dominated accretion flows (ADAFs) have been derived originally (Narayan and Yi 
1995) only for accretion onto BHs. It has been demonstrated that observed reversal events are 
reproduced by uj a ~ 0.1 — 0.4 when the transition occurs at M = M c ~ W 16 g/s. This supports 
the existence of the plasma roughly with equipartition in accretion flows around LMXBs when 
M < M c . 

The critical accretion rate M c ~ 0.la 2 MEdd where M-Edd is the Eddington accretion rate with 
10% efficiency (Narayan h Yi 1995). For a = 0.3, M c ~ 10 16 g/s is expected. This rate could 
be somewhat lower than the critical rate for the ADAF black hole systems, since cooling of the 
accretion flow is more efficient in neutron star systems due to the soft photons from the neutron 
star surface (Narayan & Yi 1995). The appearance of the quasi-spherical flow is accompanied by 
slight luminosity changes (Yi et al. 1997 and references therein) but the X-ray spectral changes 
appear to be significant (Vaughan & Kitamoto 1997). 

When the accretion flow rotation becomes sub-Keplerian, the corotation radius is shifted to a 
new location, R c = uj a R c (Yi et al. 1997). The magnetospheric radius R Q is then shifted to the 
radius R' Q which is determined by (cf. eq. 2-1) 

(R'JKf [l " (K/K) V2 } 1 = 2N c /u a N' Q (3-1) 

where N' a = u> a M(GM^R' ) 1 / 2 . The torque exerted on the star is modified as (Yi et al. 1997) 

N> = (7N>/6) [l - (8/7) (R'JKf 12 ] I [l - [K/Kf 2 ] ■ (3-2) 

In the spin-up regime, the spin-up torque is given by 

N' w \2uj a {GM,f/ 7 Rl /7 ^ /a) 1 ' 1 B 2 J 7 M % /\ (3-3) 

and the spin-up line is given by 

P* = [49m S ]^_ l7 3 / 7 a:? /7 < 8 7 Mr 6 3/7 (3-4) 

where u a ,—i = w a /0.1. The equilibrium spin period for given M and becomes longer by a 
factor uj^ 1 > 1. P* ^ 5ms would be reached for 

B,<[3x 10 8 GVl /6 (<*/7) 1/2 Mi 6 . (3-5) 

For uj a = 0.2 and M = 10 16 g/s, P* ~ [15ms]P^g 7 , which suggests that a typical MSP with 
P* = 5ms would require B» ^ 3 x 10 7 G. Any LMXBs with M < 10 15 g/s cannot produce the 
observed MSPs with B»>5x 10 7 G. 

It is usually thought that neutron stars are born with larger magnetic fields and that such 
low magnetic fields as required (P* ^ 3 x 10 7 G) for spinup of low luminosity LMXBs to MSPs 
must therefore result from field decay. Although the exact mechanism for the field decay still 
remains uncertain, in a simple heuristic model, the field decays as a result of mass accretion and 
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the strength of the magnetic field is determined by the amount of mass accreted (e.g. Taam & van 
den Heuvel 1986, Romani 1990). In this model, B* oc (MAf/AM^)" 1 , where At is duration of 
accretion and AM^ is the characteristic mass scale for field decay. Assuming that neutron stars are 
typically produced with field strengths <; 10 12 G, the spin-up of atoll sources with M ~ 10 17 g/s to 
MSPs requires At ~ 5 x 10 7 yr and B*~3x 10 8 G (and ~ 6 x 10 8 yr for the sub-critical accretion 
case with M ~ 10 16 g/s and B t <3x 10 7 G). Thus for the atoll source case (with M ~ 10 17 g/s), 
the field decay model requires that the characteristic mass AM^ ~ 3 x 10~ 5 M Q for x ~ 1. 

However, there is now evidence against field decay (in LMXBs) to values below ~10 8 G: the 
newly discovered LMXB-MSP SAX J 1808.4-3658 (Wijnands and van der Klis 1998) directly shows 
this LMXB contains a NS with a likely magnetic field B ~ (2 — 14) x 10 s G (although these same 
authors speculate that perhaps the lack of detection of pulsations in other LMXBs suggests they 
may indeed have lower B fields). Additional evidence against field decay (to below 10 8 G) in even 
very old LMXB systems comes from the detection of kHz QPOs and possible pulsations (at the 
QPO difference frequency of 275±8 Hz) in the relatively luminous atoll source 4U1820-30 in the 
globular cluster NGC 6624 (Smale et al. 1997, Zhang et al 1998). Since the NS in this system is 
a Pop II NS, and therefore likely to have an age ~10 10 yr, and yet the kHz QPOs (and possible 
pulsations) in this object with M ~ 10 17 <7/s suggest a magnetic field B ^ 3 x 10 8 G, it appears 
that even continued accretion has not reduced the B field below this value. 

Thus the required B field values of MSPs actually observed, the required spin-up timescales, 
and the lack of evidence that field decay continues below ~10 8 G combine to effectively rule out 
any persistent low luminosity sources as a possible source of MSPs. In fact, spin-up of atoll sources 
to MSPs also becomes questionable. This rather radical conclusion could be avoided only when 
a/7 > 1 or when M c <^ 10 16 g/s. The former possibility is not likely (Wang 1995, Aly Sz Kuijpers 
1990) and the latter is not supported by the recent work (e.g. Narayan & Yi 1995, Yi et al. 1997, 
and references therein). 

In Figure lb, we show the spin-up lines assuming that the accretion flow becomes sub- 
Keplerian below the assumed critical accretion rate 3.5 x 10 16 g/s with uj a = 0.2 (Yi et al. 1997). 
This figure confirms our conclusion that even at a relatively high accretion rate ~ 10 16 <7/s, 
MSPs would be produced only when B,^3x 10 7 G. For such low B*s, most of LMXBs with 
M < 6 x 10 15 g/s would not emerge out of the so-called pulsar graveyard below the death line 
(e.g. Ritchings 1976) during spin-up. They would be detected as MSPs only after arriving at 
periods P, ^ 3 x 10~ 2 s (i.e. solid lines above the the death line). It is interesting to point out 
that in the vs. B* plot, there could be a gap due to the accretion flow transition as implied 
by the spin-up line for the assumed M c ~ 3.5 x 10 16 <7/s. For M > M c the Keplerian spin-up line 
in Figure la applies whereas for M < M c (Figure lb) the sub-Keplerian spin-up line would be 
shifted discontinuously downward, leaving a region between the two lines as a gap. If there is a 
significant population of low luminosity LMXBs spun-up with the sub-Keplerian rotation, they 
could produce MSPs above the death line typically with P* ~ 5ms and -B* ~ 10 7 C These MSPs 
would populate the region marked by the solid box in the bottom left corner in Figure lb. Their 
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intrinsic electromagnetic power is expected to be only ~ B%R®Q^/6c? ^ 2 x 10 30 erg/s, which is 
roughly two orders of magnitude lower than that of the typical MSPs and is therefore likely well 
below the ~lmJy-kpc 2 radio luminosity limit of MSPs detected by Lyne et al (1997). It is not 
clear whether a significant number of such MSPs could be found in reality given the longer spin-up 
time scales (<; 10 9 yr) for lower M's. 

4. MSPs from Reccurent Soft X-ray Transients 

Some LMXBs show low luminosities (i.e. M < M c ) between recurrent outbursts (Verbunt et 
al. 1994) or during some periods as persistent sources (e.g. Tanaka & Shibazaki 1996). During 
high M, the spin-up time scale and the equilibrium spin period could become considerably shorter 
than those of the low M period. If the high M and low M periods alternate, the LMXB systems 
would be constantly away from their equilibrium spin periods. The typical spin-up time scale 
t spin-up = P*/P* depends on M as 

t spi n-up = L/P*N « 2(M,Rl)^ 7 /( 1 /a) 1 / 7 (GM 2 f/ 7 B 2 J 7 P, 

« [5 x itfyr^/a/l^M^B^p-l (4-1) 

where P* : \ = P^/Wms. This time scale applies to both the Keplerian accretion disk and the 
sub-Keplerian flow. If recurrent outbursts have M as high as ~ 10 18 <7/s, the short equilibrium 
spin period corresponding to this high M would be reached if the outburst phase cumulatively 
lasts for <; 10 7 yr while accreting a mass of ~ 0.1M & . Since the low M < W 15 g/s phase would 
have a long spin-up time scale > 10 9 yr, the recurrent outbursts would determine the eventual 
spin period if the duration of the outburst is & 10~ 3 of the quiescent duration. For instance, for 
an outburst recurrence time scale ~ lOyr, the required duration of outbursts is ~ lmonth, which 
is typical of the "soft x-ray transients" (SXTs) containing black hole primaries (e.g. Tanaka & 
Shibazaki 1996). 

SXT systems containing neutron star primaries (NS-SXTs) could thus in principle solve the 
MSP formation problem. The plausibility of this channel of MSP formation critically depends 
on actual detections of NS-SXTs in sufficient numbers. Presently, there is little evidence which 
indicates the existence of a large population of such transient sources, and in fact van Paradijs 
(1996) has shown that among the transients BH-SXTs are likely to be favored over NS-SXTs. Of 
the known NS-SXTs, Aql X-l has been found in recent RXTE observations (e.g. White & Zhang 
1997) to have a ~ 500i^2! pulsation, which could put Aql X-l as a promising MSP progenitor. 
However, since it has a ~ lyr recurrence time scale and ~1 month outburst duration and so 
is more similar to the quasi-persistent source 4U1608-52 and the atoll sources, for which the 
birthrate problem exists, it may not represent the required class of NS-SXT progenitor. In fact 
only one bona-fide NS-SXT system is known: Cen X-4, which shows deep quiescence between 
strong outbursts every ~ 10-20 years. More detections of NS-SXTs (as discussed below) are 
definitely required for them to be regarded as likely MSP progenitors. 
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Recently, Gotthelf and Kulkarni (1997) have reported the detection of an unusual burst 
source in the globular cluster M28. They attributed the unusually low burst luminosity to the 
low mass accretion rate onto a magnetized neutron star. The possible type I burst could occur 
in the magnetized accretion column which naturally increases the effective mass accretion rate 
and density on the stellar surface. Even if such sources eventually show up in great numbers as a 
potentially important new class of NS-SXTs, the deduced magnetic field strength (based on the 
non-detection of pulsations and hence extremely slow spin) is much too large to be compatible 
with the typical field strength needed for MSP progenitors. 

5. Summary and Discussion 

We have found that low luminosity LMXBs are hard to spin-up to MSPs unless they have 
magnetic fields substantially lower 5 x 10 7 G) than those of the observed MSPs with B ^ 10 8 
G. Thus the LMXB phase would have to be followed by a phase of field amplification to produce 
the currently observed MSPs. Even when the magnetic fields are sufficiently weak, the long 
spin- up time scale ~ 10 10 yr for low mass accretion rates <C W 15 g/s poses another problem. If, as 
is likely for ADAF models (Narayan & Yi 1995), the accretion flow becomes quasi-spherical and 
sub-Keplerian for M < g/s, it could rule out any low luminosity population and potentially 
a significant fraction of atoll sources as a possible source of the MSPs. Thus we conclude that 
the "standard model" for the formation of MSPs from LMXBs requires that they form from a 
larger population of normally quiescent SXT systems containing NSs with B ^ 10 8 G. If these 
sources have fractional duty cycle f = 0.01 (i.e. accrete at ~ 10 18 g/s for ~1 month every ~10 
yr), they could provide the observed MSP population if their parent population is ~ 10 x the 
observed persistent atoll source population (since, although only ~ 10 7 yr is needed for spin-up, 
the ~ 1% duty cycle means that an effective lifetime needed is 10 9 years, or the same as that for 
the persistent sources). More generally, a population of quiescent NS-SXTs with peak accretion 
rates Mis ( m units of 10 18 g/s) and on-time duty cycle f_2 (in units of 0.01) could only solve the 
MSP birthrate problem if at any given time there are Nns-sxt ~ 10/(Mi§f_2) yr -1 in outburst. 
Thus there would need be ~10 NS-SXTs like Cen X-4 (with M 18 ~ 1 and f_ 2 ~ 1) or ~100 like 
Aql X-l (with Myj ~ 1 and f_2 ~ 10) in their "on" state at any given time, which is clearly not 
consistent with observations. Although the discovery of SAXJ1808. 4-3658 may point to a larger 
population of Aql X-l type NS-SXTs, the apparent rate of these new bursters discovered by RXTE 
or BeppoSAX thus far appears much too low: perhaps at most ~3 "new" (recurrent) bursters 
that are "on" at any given time. A recent example is the RXTE re-discovery of the burster 
XTEJ1806-246 (Marshall and Strohmayer 1998), which is very likely the burster 2S1803-245 
originally detected with SAS-3. Thus the required large population of NS-SXTs are simply not 
seen (it would be comparable in number to the known number of persistent atoll sources), and so 
we conclude there is as yet no compelling evidence for MSP formation from NS-SXTs. 

A new class of MSPs could possibly form from an as yet undetected population of persistent, 
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low luminosity LMXBs. These would produce a low magnetic field (< 10 7 G), radio weak MSPs 
with typical radio luminosities ~ W 2S erg/s if the radio efficiency is ~ 1%. If the magnetic field 
is much weaker than ~ 10 7 G, the MSPs would reside in the region below the pulsar death line, 
which makes the detection of these pulsars virtually impossible. If a significant number of low 
luminosity LMXBs spin-up to MSPs, they could appear as a group of radio-weak MSPs distinct 
from the known MSPs. 

Given the difficulties of the "standard model", even for NS-SXTs, we consider finally the 
accretion-induced collapse (AIC) model for "direct" production of MSPs without the LMXB 
progenitor phase (Michel 1987, Grindlay & Bailyn 1988, Bailyn & Grindlay 1990). In the 
AIC scenario, white dwarfs collapse to neutron stars after reaching the critical Chandrasekhar 
limit ~ 1.4Mq through accretion. It is believed that AIC occurs for high mass accretion rates 
> 2.5 x W 18 g/s (e.g. Nomoto & Kondo 1991, Livio k Truran 1992). The pre-AIC magnetic field 
on the WD will play an important (and poorly determined) role in the formation of MSPs. For 
magnetic WDs (B ?Z 10 6 gauss) and initially slowly rotating white dwarfs accreting at a rate in 
the range ~ 1 x 10 18 — 7 x W ls g/s and assuming that the magnetic fields are flux-frozen and 
amplified by compression during collapse, Yi & Blackman (1997) have found that the magnetic 

a It: 

field of AlC-produced pulsars would be correlated with the pulsar spin as fi* 4 ~ rjB^ 1 | (cf. 
Narayan & Popham 1989) where n* )4 = Ot/lO 4 ^ 1 , £*,n = 5*/10 n G, and 7] is a constant of 
order unity. In Figure lb, the dotted line is the hypothetical birth line of strongly magnetized 
pulsars from AIC. This correlation indicates that the field strength of pulsars with ~ 10 3 s _1 
would be typically close to ~ 10 13 G, which are many orders of magnitude larger than the observed 
MSPs. Unless these fields decay rapidly, or if magnetic flux is not conserved during collapse in 
AIC (e.g. if instead field is annihilated), the AlC-produced MSPs would be very bright and lose 
their rotational kinetic energy at a rate L*,_bm ~ 10 44 f2^4 erg/s with the characteristic emission 
time scale ~ I*Q%/2L* t EM ~ 2 x 10 8 f^4,s. Such a short time scale indicates that the pulsars 
could evolve almost instantaneously to long P*'s. They could be an abundant source of LMXBs if 
accretion resumes after ~ 10yr of intense electromagnetic activity. These AlC-produced LMXB's 
would produce MSPs similar to the observed ones only when the magnetic field somehow decays 
by several orders of magnitude during the LMXB phase. However, by requiring a subsequent 
LMXB phase, they obviously do not solve the observed LMXB vs. MSP birthrate problem and so 
we do not consider this further (for MSP production). 

For the case of weakly magnetized (^ 10 4 ~ 5 G) white dwarfs, AIC could directly produce 
low-field MSPs but only if the collapse itself is not halted by rotation. To produce the observed ~2 
ms spin periods, the pre-AIC WD must have spin period ~20 s (assuming conservation of angular 
momentum in the collapse from a 1.4M WD of radius ~ 10 3 km to a 10 km radius NS). However, 
for low magnetic field, these stars are not braked by the magnetic fields and so the spin period 
of the pre-collapse WD could instead be spun up to the minimum ((47r 2 -R 3 ^/GM^) 1 / 2 ~1s) 
period where collapse would be halted unless the excess angular momentum is first removed. 
Thus for AIC to be completed without producing the so-called rotationally supported "fizzlers" , 
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angular momentum must be efficiently removed during or before collapse (e.g. Tohlin 1984). This 
implies that AIC would produce the MSPs only if white dwarfs stop spinning-up despite continued 
accretion. Thus the possibility of spin-down near the break-up period, despite continuous accretion 
(Popham & Narayan 1991, Paczynski 1991, Bisnovatyi-Kogan 1993), is an intriguing possible 
solution to the MSP birthrate problem. 

We conclude that if an AIC process produces the bulk of the observed MSPs, it is necessary 
to shed angular momentum (low B field WD case) or magnetic field (high B field WD case) either 
prior to or during the collapse. Since only the former process directly produces a MSP without an 
LMXB, it is preferred. However both processes may occur and in fact may relate the problems of 
MSP formation with that of the origin of gamma-ray bursts, as we discuss in a forthcoming paper 
(work in preparation). 

We thank the referee, S. Kulkarni, and J. McClintock for useful comments on the manuscript. 
This work was supported in part by NASA/LTSA grant NAG5-3256 (to JEG) and by the SUAM 
Foundation (to IY). 
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Figure Captions 

Figure la: Spin-up lines in the vs. P* plane. The x's are the observed pulsars adopted from 
Taylor et al. (1995). The dashed line is the pulsar death line, B* ~2x 10 U P^. The solid lines are 
the spin-up lines for M, W 18 g/s, M c = 3.5 x W w g/s, W 16 g/s, 6.3 x W 15 g/s, and W 13 g/s from 
top to bottom, respectively. The thick solid lines are relevant for typical atoll sources. 

Figure lb: Similar to Figure la with the sub-Keplerian rotation (uj a = 0.2) for M < M c . The 
spin-up line for a given M < M c moves down due to the sub-Keplerian rotation. iFvom top to 
bottom, M's are 3.5 x W 16 g/s, W w g/s, 6.3 x 10 15 g/s, and 10 13 g/s, respectively. The solid box 
in the bottom left indicates the location of the "typical" MSPs spun- up by the sub-Keplerian 
accretion flows. The dotted line is a possible birth line from AIC for M = 3.2 x 10 18 g/s. The 
thick solid lines are relevant for typical atoll sources. 



